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1 Introduction 

In aerospace applications, magnetic suspension systems may be required to operate over large variations 
in air-gap. Thus the nonlinearitics inherent in most types of suspensions will have a significant effect. 
Specifically, large variations in operating point may make it difficult to design a linear controller which 
gives satisfactory stability and performance over a large range of operating points. 

One way to address this problem is through the use of nonlinear compensation t echniques such as 
feedback linearization. In the last, decade, this area has seen a great deal of research eflort among system 
theorists [()], [7], [8], and has been successfully used in flight control applications [9], but has yet to be 
widely used in the area of magnetic suspensions. 

Nonlinear compensators have received limited attention in the magnetic suspension literature. In 
[11] the control system design for suspension of a 1-ton, 4-passenger vehicle is presented. A flux sensor 
is located on the pole-face of the suspension magnet. A minor feedback loop is closed on flux which 
linearizes the dependence of magnetic force on position. The force goes as the square of flux, so an 
additional square root linearization is required. The idea of using an analog multiplier to compute the 
ratio i/r. is mentioned, which would allow elimination of the flux sensor. However, this technique is 
reported as prone to drift and noise, and was abandonded. This flux-feedback scheme is described in 
more detail in [12]. In [13] a nonlinear correction law is used to correct tin 1 inverse square law magret 
behavior in a flywheel suspension. The nonlinear compensation is implemented with analog multiplier 
and square root circuits. Microprocessor- based linearizing transformations are reported in [14] in she 
context of a demonstration system. 

In recent, years, progress has been made in the theory of nonlinear control systems, and in the sub- 
area of feedback linearization. Here, [6] is of fundamental importance in that it presents the conditions 
under which a system may be linearized. In [7] these results are globally extended, and in [8] the theory 
is developed for rnulti-input/multi-oiitpiit systems. In a subsequent section the results of [6] are applied 
to a third order suspension example. However, for a simple system, it is often possible to construct 
the linearizing transformations by inspection. We start then by demonstrating the idea of feedback 
linearization using a second order suspension system. In the context of the second order suspension, 
sampling rate issues in the implementation of feedback linearization are examined through simulation. 
The system which is studied is presented in the next section. 


2 Nonlinear Suspension Model 

In this section, the open-loop dynamics for a simple oue-degreo-of-freedom suspension ere presented. ' his 
system exhibits the essential issues faced in the design of tractive type suspensions, that is, suspeie >ns 
which operate as variable reluctance devices. The example is drawn from [10] pgs. 22-23, 84-80, ind 
193-200. The only change is that the system is inverted such that gravity acts to open tin* air-gap. I his 
system is shown in Figure 1. 

The details of the electromagnetics are worked out in [10]; for the present purposes, the import int 
details are the coil voltage 
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Figure 1: Single degree of freedom suspension. 



and the iom* on the plunger 
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where the first term is the electromagnet force, the second term is the gravitational force on the plunger, 
and the third term is a disturbance force acting on the plunger in the direction the electromagnet force. 
The nonmagnetic sleeve is assumed to exert no frictional forces on the plunger. 

If we define 
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then the state equations for the open-loop suspension arc 
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2.1 Second order system 

If the coil current is assumed to be the cont rol input, then the suspension slate equations are reduced 
to second order. 
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These equations will adequately model the system if the coil current is controlled by a high-bandwidth 
current loop with sufficiently high voltage-drive capabilities. In applications, it is most typical to drive the 
coil with such a current loop, as this essentially eliminates the dependence of position-loop performance 
upon the electromagnet coil resistance and inductance. 


3 Linearization of second-order suspension 

The basic idea of feedback linearization is to define transformations on the states and input(s) such 
that the nonlinear system appears linear and operating-point invariant in terms of the transformed 
representation. Then a controller can be designed for the transformed variables. This allows the closed- 
loop system stability to be made independent of operating point. 

For the second-order equations (4), a transformation on the input is all that is required to linearize 
the system. This transformation may be derived by inspection without using any formal mathematical 
machinery; this is the approach taken in [11], [12], [13], and [14]. That is, if the coil current i is made to 
vary as 


* = x i 



(■ r >) 


then the suspension is globally linearized in terms of the new input v. 'The notation for the auxiliary 
input v lias been chosen to match the notation in [b], 

Specifically, substitut ing from (5) into (d), the system slate equations become 
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These equations are linear, with an input t>, and disturbance terms g and fd- 

Here v is a signal internal to the compensator which may be thought of as a setpoint for accelerat.on 
in the direction of increasing airgap. In operation, the signal v will be computed within the compensator, 
and constrained to be less than or equal to zero. Since the magnet can only supply accelerations in he 
direction of decreasing air gap it would not be physically meaningful to ask for acceleration in the 
direction of increasing air gap by setting v greater than zero. Thus the term -v m (5) will always be 

greater than or equal to zero, and the square root will yield a real number. 

The plant appears linear in terms of the new input v. This compensation of the nonlinear term oes 
not however stabilize the plant. To stabilize the system, the nonlinear compensator is preceded by a 
linear compensator. The resulting closed-loop system is shown in Figure 2. The compensator may be 
thought of as having two parts, a nonlinear compensation section and a linear compensation section. It is 
the function of the nonlinear section to implement (5) in order to adjust i as a function of *i and v such 
that the acceleration of the plunger is equal to v. It is the function of the linear section to specify the 
value of v as a function of the error between the position setpoint and the measured position such tha 
the linearized plant is robustly stabilized and has good disturbance rejection and settling time properties. 
The signal v forms the connection between the linear and nonlinear sections of the compensator. 

This combination of linear and nonlinear compensation sections stabilizes the plant such that the loop 
dynamics are independent of operating point. Such operating point independence is the main advantage 
of using a nonlinear compensator. Note that as viewed from the input to the nonlinear section, the 
incremental relationship between v and *, is equal to 1/s 2 , independent of operating point. Thus, the 
linear compensator can be designed to control a double integrator via standard linear techniques. If l is 
desirable to reject static disturbance forces with no position error, then the linear compensator can be 
designed to include an integral term. This integral term will adjust the value of v to balance gravity an 

any low-frequency components of the disturbance /<*• i . 

In applications where large excursions or disturbance forces are anticipated, the additional complexity 
of the nonlinear compensation approach is justified. The major caveat is that we are assuming t lat t c 
suspension model is accurate. For the electromagnetics an accurate model can readily be developed, and 
thus nonlinear compensation techniques are applicable. The nonlinear compensation technique was use 
in the construction of a class demonstration system which is described below. 


4 Classroom demonstration implementing linearization 

In the Spring of 1988, the author constructed a single degree of freedom levitation system for use as 
a classroom demonstration which implemented the nonlinear compensation technique described in the 
previous section. As developed there, if the plant state equations are given by (4), then app ying t le 
nonlinear compensation law (5) results in a system which appears to be linear in terms of the intermediate 
signal v. The demonstration system uses a high-bandwidth current-drive to regulate the electromagnet 
current, and thus (4) is applicable. 

In the demonstration system, a one inch steel ball bearing is suspended below an electromagne 
consisting of 3100 turns of #22 magnet wire wound on an 1 inch diameter by 4 inch length steel core. 
The coil current is controlled by a Bose-type switching regulator, with a half-scale current switching 
frequency of 10 kHz, and a full scale current of 2 Amperes. The operating point current ls about 0T 
Amperes at a typical operating point air gap of 1 cm. The system is digitally controlled by an 8088/8087- 
based single-board computer and data acquisition system at a 400 Hz sampling rate. The contro law 
for the nonlinear compensation section uses (5) to linearize the magnetic force relationship. Ih.s allows 
the stability of the closed-loop system to be essentially independent of the operating point. The contro 
law for the linear compensation section is then developed via classical techniques applied in the discrete- 
time domain. The position of the ball is sensed optically, and nonlinearities in the sensor output versus 

position are compensated for in software. . , 

In order to apply the nonlinear compensation technique, an accurate model of the plant is required. 
For the classroom demonstration, this model is developed by measuring the force on the ball as a function 
of current and position. This measurement is accomplished by using a balance beam for measuring e 
magnetic force on the ball. A 1 inch ball bearing is glued into one end of an aluminum balance b^rn 
of rectangular tubular cross-section measuring 1 inch wide by 0.75 inches deep by 12 inches long. The 
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figure 2: Nonlinear compensation of second-order magnetic suspension system. 




beam pivots at the center on a thin wire which is held by fixed side supports. B /u " ! ^ h ese 

£ 2 : r - je : r *. u- 
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“d" f i JH that i, interacts with the optical sensor in the ,„nc fasMon . a freely »,p»<M 

bal The force relationship (2) was well lit by the experimental data with the parameter. C = *•« » 
10-< NmW »d *, = 0 Wen. The man. of the hall la 67 grama. Th« f parameters are used in the 
nonlinear compensation law (5). The only deviation from the relationship (2) is at high corren 

where the effects of magnetic saturation are apparent. . 

The optical position^ sensor is constructed a. follows. A 24 volt, 5 watt incandc*,» lamp » “ 

the source and a piece of cadmium sulfide photo-cell is used as the sensor, in what is a standard po 
sensor for magnetic suspensions. Using the balance beam described above, the sensor output « me^red 
for a number of ball positions. When the shadow-line cast on the sensor is in the central g 
sensor the sensor output is essentially linear with ball position. However, as the shadow-line approache 
the upper or lower edge of the photocell, the sensor sensitivity begins to decrease. This nonlinearity 
the relation between ball position and sensor output is corrected in software in the section of code wh. 
Lpu^hTJLr voltage^ The corrected position measurement then linear wit £«., hall position. 
It is this corrected position measurement which is passed to the rest of the control loop. 

The position sensor was found to have several defects which limit the system performance First, 
the Incandescent bulb output decreases significantly as a function of time. This is be keyed to be due 
o the evapo I ion of the filament. Material driven off of the filament ,s deposited on the inside of he 
rf^ enTpe hereby decreasing the bulb brightness. The second problem is that the -dm.um sulfide 
sensor is sensitive to any light falling on its surface, independent of the source. Thus ambient g ■ g - 

indistinguishable from the light emitted by the bulb. correction 

Both of these effects cause problems in the nonlinear compensation law (5) . 

of the til ZZritim. First, the decrease in bulb intensity and any changes m average ambient 
lieht act as offset terms which drive the system to incorrect points on the sensor correc ion curve 
in the magnet nonlinearity correction law (5). This offset deteriorates the system stab. 
the ambient light has a large component at twice the power line frequency, especia > 
fluorescent lighting. This signal at 120 Hz acts as a large noise source which causes error motions m 

ba ^The S above problems can be solved as follows. First, the light source needs to be made more constant 
with time This can be achieved by using a more specialized incandescent bulb, or y swi e g 
semiconductor light source such as an infra-red light emitting diode The ambient lighting offset and 
noise problems can be solved by either or both of two approaches which are classical, rhe firs is to make 
the system narrow-band. Commonly available 1R diodes emit a relatively narrow-band optical signal 
^r dtd" are narrower. In this cal, an optica, band-pass filter can be placed in front of he sensor so 
that only the emitted frequencies are sensed, and the ambient lighting is greatly attenuated The second 
approach is to switch the light-source on and off at a high frequency and use synchronous detection to 
reject signals which are not at the same frequency and phase as the source. The requency o swi c mi 
must be*inade much higher than the cross-over frequency of the position control loop perhaps on the 
order 5 .0 kHz switching frequency. This rate is easily within the capabilities of av^We e ec jomc. 
The results derived in the previous section for the nonlinear compensation laws assume that these are 
uL. For discrete-time implementation, the issue ot sampling rate becomes 

important. This issue is investigated in the next section. 
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4.1 Sampling rate issues 

Due to the complexity of the transformations it is most likely that a linearizing compensator will be 
implemented in discrete time. As an introduction to one issue involved in discrete-time implementation, 
the effect of sampling rate on the second-order suspension system (4) is investigated by simulation. For 
this example, the suspension parameters have been been given the values developed for the class demon- 
stration system described above. These values are M = 67 grams, and C = 4.43 x 10“ 4 Nm 2 /Amp 2 . 

The system is simulated assuming a nonlinear compensation law of the form (5). The four graphs 
shown in Figure 3 indicate the system behavior when a net 0.05g acceleration (v = —0.05 in (5) ) is 
specified. The lines labelled ‘ideal' show that if the nonlinear compensation was perfectly implemented, 
the force on the ball would be constant, and the graph of velocity vs. time would be a straight line. 
However, with any finite sampling rate this is not the case. The system is open-loop unstable, and 
uncontrolled between sampling instants. Thus it ‘runs away' during the interval in which the control 
current is held constant. The graphs show the result of this process for sampling rates of 1 kHz and 200 
Hz. To get reasonable behavior, it can be seen that a sampling rate on the order of 1 kHz is required. 
In the class demo, due to computational speed limitations a 400 Hz sampling rate is used. This is found 
to be adequate as long as the ball is not allowed to approach too close to the pole face. 

Another way to look at the effect of sampling rate is to examine the system behavior under closed-loop 
position control. To this end, a linear proportional plus lead compensator is designed in discrete-time to 
stabilize the nominal plant which would result if the nonlinear compensation were perfect. That is, in the 
ideal case, the nonlinear compensated system appears as a double integrator independent of operating 
point. In the finite sampling time implementation, the quality of this approximation deteriorates as the 
air gap closes. This can be seen in Figure 4 which displays simulated step responses for the closed-loop 
system at four nominal operting points and for the two sampling rates. Note that the system with 200 Hz 
sampling goes unstable at the 0.5 cm and smaller air gaps, whereas the behavior of the 1000 Hz sampled 
system only begins to deteriorate when the air gap approachs 0.3 cm. The unstable response for 200 Hz 
sampling is not shown for the 0.3 cm air gap. 

This example shows that the practical implementation of linearizing transformations may require 
very high sampling rates. Also, what may be considered a satisfactory sampling rate depends on the 
range of operating points which are encountered in system operation. Certainly, the issue of discrete-time 
implementation merits further study. 

Experience with this simple nonlinear compensation system provided the impetus toward an under- 
standing of feedback linearization techniques in more generality. A description of the application of 
feedback linearization to the third order suspension system is given in the next section. 


4.2 Linearization of third-order suspension 

For more complex plants it may be difficult to develop linearizing transformations by inspection. The 
results of [6] provide a general approach to this problem. A good introduction to these ideas is presented 
in [4] and [5]. These references assume no more than an undergraduate background in control theory 
and are thus a good place to start for someone new to this area. 

Without reviewing the results from the above references, if the plant satisfies a controllability con- 
dition and a condition on the existence of solutions to a set of partial differential equations, then trans- 
formations z\ = Ti(x), . . . , z n = T n (x),v = T n + 1 (ar,u) can be constructed such that in the z-v space 
the system appears linear. Here, x is the state- vector of the nonlinear system, 2 is the state-vector of 
the linearized system, and n is the system order. Under these transformations, the nonlinear system is 
mapped to the controllability canonical form 
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For the system (3), the required conditions are satisfied, and the results of [6] yield the linearizing 
transformations Z\ — X\, = ^ 2 » z 3 — — ) 2 , and Z 4 = ~ u). Thus the system 
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Fimce 3 The open-loop system with nonlinear compensation showing its performance with sampling 
S t If 1 and 5 milliseconds. Position . is in cm separation from the pole face, velocrt, . in 
cm/sec, force on the ball is in dynes, and current I is in amps. 




1: igure 4: 1 he closed-loop system with nonlinear compensation showing its performance with sampling 
periods h of 1 and 5 milliseconds as the operating point position is moved towards the pole face. Position 
x is in cm separation from the pole face. The nominal operating points are 2, 1, 0.5, and 0.3 cm 
respectively. 


the suspended member. Thus it makes physical sense that the suspension w,| l appear lnearmj3 ; Th 
implementation, the voltage drive u must be computed m terms of t>: 

( 8 ) 


u = _**!Z}l + iR. 


Since « drives the derivative of 23 , we can think of t> as being a setpoint for the slope of the acceleration. 



Direct substitution will verify that the transformations Ti - x,, T 2 - ^x,x 2 , 


T 3 = 2 x 2 
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and 




(10) 


position x x and current *. Both of these quantities may be readily measured. 


5 Conclusions 


As we have seen in the magnetic suspension examples, the technique of feedback linearization w of 

errors. Thi. is an »» which has .bo b e en a top,, of current research 151 In o . ‘ 

compensation laws for this and other multivariable suspension systems. 


6 Acknowledgements 

This work forms part of a thesis submitted by the author to the Department of Electrical Engineer mg 
and Computer Science .1 MAT. in partial fulfillment of the requirements for the °f j^ctc r rf 

ES^MWSS c “ r T m.: m 

Mitiguy, and Hard, Prince at MAT. Lincoln Laboratory Thi. work was .pon.ored by the Department 
of the Air Force, and performed at the M.l.T. Lincoln Laboratory. 


7 References 

111 Slocum A If and Eisenh.ure, D.B., ‘Design Consideration. for Fllra-Precision Magnetic Bearing 
Supported Slide." NASA Magnetic Su.pen.ion Technology Conference, Hampton, Va. Feb. 2-4, 


[2] lYumper, D.L., “Magnetic Suspension Techniques for Precision Motion Control,” Ph.D. Thesis, Dept, 
of Elec. Eng. and Comp. Sci., M.I.T., Camb., Mass., Sept., 1990. 

[3] Trumper, D.L., “Five- Degree-of- Freedom Control of an Ultra- Precision Magnetically-Suspended Lin- 
ear Bearing,” NASA Workshop on Aerospace Applications of Magnetic Suspension Technology, Hampton, 
VA, Sept. 25-27, 1990. 

[4] Slotine, and Li, preliminary notes for text on nonlinear control, to be published Jan., 1991. Used as 
course notes for MIT class 2.152. 

[5] Spong, M., and Vidyasagar, M., Robot Dynamics and Control , Wiley, 1989, pp. 259-283. 

[6] Su, II., “On the linear equivalents of nonlinear systems,” Systems and Control Letters, vol. 2, No. 1, 
pp. 48-52, July, 1982. 

[7] Hunt, L.R., Su, R., and Meyer, G., “Global Transformations of Nonlinear Systems,” IEEE Transac- 
tions on Automatic Control, Vol. AC-28, No. 1, Jan., 1983. 

[8] Hunt, L.R., Su, R., and Meyer, G., “Design for Multi-Input Nonlinear Systems,” in Differential 
Gcometf'ic Control Theory , Brockett, Millman, and Sussmann, eds., Progress in Mathematics, Vol. 27, 
pp. 268-298, Birkhauser, Boston, 1983. 

[9] Meyer, G., Su, R., and Hunt, L.R., “Application of Nonlinear 'lYansformations to Automatic Flight 
Control,” Automat ica, Vol. 20, No. 1, pp. 103-107, 1984. 

[10] Woodson, H.H., and Melcher, J.R., Elcctivmechanical Dynamics - Tart 7, John Wiley and Sons, 

1968. 

[11] Jayawant, B V., Hodkinson, R.L., Wheeler, A.R., and Whorlow, R.J., “Transducers and Their 
Influence in the Design of Magnetically Suspended Vehicles,” I.E.E. Conf. on Control Aspects of New 
Forms of Guided Land Transport, I.E.E. Publication No. 117, Aug., 1974, pp. 200-206. 

[12] Jayawant, B.V., Sinha, P.K., Wheeler, A.R., and Whorlow, R.J., “Development of 1-ton Magnet- 
ically Suspended Vehicle Using Controlled D.C. Electromagnets,” Proc. I.E.E., Vol. 123, No. 9, Sept 
1976. 

[13] Groom, N.J., and Waldeck, G.C., “Magnetic Suspension System for a Laboratory Model Annular 
Momentum Control Device,” AIAA Paper No. 79-1755, 1979. 

[14] Traxler, A., Meyer, F., and Murbach, IF, “Fast Digital Control of a Magnetic Bearing with a 
Microprocessor,” International Kongress Mikroclektronik, Munich, Nov 13-15, 1984, pp 94-102. 



